Using bioinformatics and systems biology methods to identify the mechanism of interaction between COVID-19 and nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is considered a risk factor for severe COVID-19, but the mechanism remains unknown. This study used bioinformatics to help define the relationship between these diseases. The GSE147507 (COVID-19), GSE126848 (NAFLD), and GSE63067 (NAFLD-2) datasets were screened using the Gene Expression Omnibus. Common differentially expressed genes were then identified using a Venn diagram. Gene ontology analysis and KEGG pathway enrichment were performed on the differentially expressed genes. A protein-protein interaction network was also constructed using the STRING platform, and key genes were identified using the Cytoscape plugin. GES63067 was selected for validation of the results. Analysis of ferroptosis gene expression during the development of the 2 diseases and prediction of their upstream miRNAs and lncRNAs. In addition, transcription factors (TFs) and miRNAs related to key genes were identified. Effective drugs that act on target genes were found in the DSigDB. The GSE147507 and GSE126848 datasets were crossed to obtain 28 co-regulated genes, 22 gene ontology terms, 3 KEGG pathways, and 10 key genes. NAFLD may affect COVID-19 progression through immune function and inflammatory signaling pathways. CYBB was predicted to be a differential ferroptosis gene associated with 2 diseases, and the CYBB-hsa-miR-196a/b-5p-TUG1 regulatory axis was identified. TF-gene interactions and TF-miRNA coregulatory network were constructed successfully. A total of 10 drugs, (such as Eckol, sulfinpyrazone, and phenylbutazone) were considered as target drugs for Patients with COVID-19 and NAFLD. This study identified key gene and defined molecular mechanisms associated with the progression of COVID-19 and NAFLD. COVID-19 and NAFLD progression may regulate ferroptosis through the CYBB-hsa-miR-196a/b-5p-TUG1 axis. This study provides additional drug options for the treatment of COVID-19 combined with NAFLD disease.

Identification of critical genes and molecular pathways in COVID-19 myocarditis and constructing gene regulatory networks by bioinformatic analysis.

BACKGROUND: There is growing evidence of a strong relationship between COVID-19 and myocarditis. However, there are few bioinformatics-based analyses of critical genes and the mechanisms related to COVID-19 Myocarditis. This study aimed to identify critical genes related to COVID-19 Myocarditis by bioinformatic methods, explore the biological mechanisms and gene regulatory networks, and probe related drugs. METHODS: The gene expression data of GSE150392 and GSE167028 were obtained from the Gene Expression Omnibus (GEO), including cardiomyocytes derived from human induced pluripotent stem cells infected with SARS-CoV-2 in vitro and GSE150392 from patients with myocarditis infected with SARS-CoV-2 and the GSE167028 gene expression dataset. Differentially expressed genes (DEGs) (adjusted P-Value <0.01 and |Log2 Fold Change| ≥2) in GSE150392 were assessed by NetworkAnalyst 3.0. Meanwhile, significant modular genes in GSE167028 were identified by weighted gene correlation network analysis (WGCNA) and overlapped with DEGs to obtain common genes. Functional enrichment analyses were performed by using the "clusterProfiler" package in the R software, and protein-protein interaction (PPI) networks were constructed on the STRING website (https://cn.string-db.org/). Critical genes were identified by the CytoHubba plugin of Cytoscape by 5 algorithms. Transcription factor-gene (TF-gene) and Transcription factor-microRibonucleic acid (TF-miRNA) coregulatory networks construction were performed by NetworkAnalyst 3.0 and displayed in Cytoscape. Finally, Drug Signatures Database (DSigDB) was used to probe drugs associated with COVID-19 Myocarditis. RESULTS: Totally 850 DEGs (including 449 up-regulated and 401 down-regulated genes) and 159 significant genes in turquoise modules were identified from GSE150392 and GSE167028, respectively. Functional enrichment analysis indicated that common genes were mainly enriched in biological processes such as cell cycle and ubiquitin-protein hydrolysis. 6 genes (CDK1, KIF20A, PBK, KIF2C, CDC20, UBE2C) were identified as critical genes. TF-gene interactions and TF-miRNA coregulatory network were constructed successfully. A total of 10 drugs, (such as Etoposide, Methotrexate, Troglitazone, etc) were considered as target drugs for COVID-19 Myocarditis. CONCLUSIONS: Through bioinformatics method analysis, this study provides a new perspective to explore the pathogenesis, gene regulatory networks and provide drug compounds as a reference for COVID-19 Myocarditis. It is worth highlighting that critical genes (CDK1, KIF20A, PBK, KIF2C, CDC20, UBE2C) may be potential biomarkers and treatment targets of COVID-19 Myocarditis for future study.